COMPUTING LEVEL ONE HECKE EIGENSYSTEMS (MOD p) 



CRAIG CITRO AND ALEXANDRU GHITZA 

Abstract. We describe an algorithm for enumerating the set of level 1 Systems of Hecke 
eigenvalues arising from modular forms (mod p) . 



1. Introduction 

One of the cornerstone results of the modern arithmetic theory of modular forms associates 
to every level 1 Hecke eigensystem mod p a unique odd semisimple 2-dimensional Galois 
representation (mod p) unramified outside p. This follows from the corresponding results 
of Deligne (and Serre, and Eichler-Shimura) for eigenforms over Z; a more direct approach 
that avoids using the füll machinery of Deligne's characteristic zero theorem can be found 
in Proposition 11.1 of [Gro90]. 

Serre's conjecture (now a theorem of Khare-Wintenberger) says that all Galois representa- 
tions described above arise from level 1 eigensystems. In Section 8 of [Kha07], Khare recalls 
the well-known fact that the set of level 1 eigensystems (mod p) is finite of cardinality 0(p 3 ) 
as p — > oo, and he outlines an argument due to Serre showing that this cardinality is o(p 2 ) 
as p — > oo. Khare adds: u It will be of interest to get quantitative refinements of this" and 
guesses that the cardinality is in fact asymptotic to p 3 /48 as p — > oo. In his PhD thesis, Cen- 
teleghe studies this question and proposes a precise conjecture for the asymptotic behavior 
of the number of representations of fixed conduetor N (see Conjecture 4.1.1 in [Cen09]). 

The present paper describes an efheient algorithm for enumerating the set of level 1 eigen- 
systems (mod p), and hence also the set of odd semisimple 2-dimensional Galois representa- 
tions (mod p) unramified outside of p. The theoretical framework underlying our approach 
is based on Tate's theory of theta cycles. We use two alternative computational methods: 
the Victor Miller basis for modular forms of level 1, and modular symbols over finite fields. 

In his recent preprint [CcnlO], Centeleghe attacks the problem of counting the number of 
irreducible Galois representations by an ingenious approach that requires Computing with 
a single Hecke Operator for each prime p. Unfortunately, this method only gives a lower 
bound on the number of representations. It is worth noting, however, that this lower bound 
is generally very close to the known upper bound, and in many cases (164 of the 299 cases 
considered in [CenlO]) allows one to deduce the exaet number. See Section 7 for more on 
the relationship between Centeleghe's work and ours. 
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2. Review of modular forms mod p 

We recall the definition of modular forms mod p of level 1 and of their Hecke Operators. 

Let Mfc(C) denote the complex vector space of holomorphic modular forms of weight k 
and level 1. There is a C- linear map that associates to each modular form its g-expansion 
at the (only) cusp oo: 

oo 

Q:M»(C)— ►C[[g]] l f — ► f(q) = ^a n g™. 

n=0 

By the g-expansion principle (Theorem 1.6.1 in [Kat73]), this map is injective. 
We dehne the Z-module of forms with integer coefficients by 

M k (Z) = Q- 1 (Z[[g]]) 

and, for any Z-module R, we dehne the .R-module of forms with i?-coefficients by 

M k {R) = M fc (Z) ®z R. 

We will dehne 1 the space of modular forms mod p of level 1 and weight k to be M k = 
M k (¥ p ). These are obtained by reducing modulo p the g-expansions of the modular forms 
with integer coefficients. 

2.1. Eisenstein series mod p. There are two normalisations for Eisenstein series in char- 
acteristic zero. The hrst makes the coefficient of q be one: 

B °° 

(1) G k = -^ + J2°k-i(n)<l n - 

n=l 

The second one makes the constant coefficient be one: 

1k 2k J^U 

(2) E k = -—G k = 1 - — r^i(n) g ". 

B k B k ^ 

We dehne Eisenstein series (mod p) by reducing the characteristic zero Eisenstein series 
modulo p. The hrst normalisation is problematic for primes dividing the denominator of 
B k /(2k); by the von Staudt-Kummer congruences (see Lemma 4 in [SD73]), this happens if 
and only if k is a multiple of p — 1. 

Convention: To simplify notation, we will always write G k to denote the Eisenstein series 
(mod p) of weight k, keeping in mind that it is the reduction modulo p of the q-expansion 
in (1) if k is not a multiple of p — 1, and the reduction modulo p of the q-expansion in (2) 
if k is a multiple of p — 1 . 

Since we will soon restrict our attention to forms of weight < p + 1, the latter Situation 
will only occur for the Hasse invariant A, which is the reduction modulo p of E p _i. The von 



^Morally, the appropriate definition of modular forms mod p is intrinsic, as global sections of line bundlcs 
over the moduli Stack of elliptic curves over ¥ p (see Section 1.1 in [Kat73], Section 10 in [Gro90], or Section 2.1 
in [Edi92]). The naive definition we use is equivalent in level 1 for p > 5, by Theorem 1.8.2 and Remark 1.8.2.2 
in [Kat73]. 
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Staudt-Kummer congruences teil us that, apart from the constant coefficient, all coefficients 
of E p _i are divisible by p, so the g-expansion of A is simply A(q) = 1 G F p [[g]]. 

2.2. Operators. The spaces Mk are equipped with a number of interesting linear maps. 
We will dehne them in the most economical way, by describing their effect on g-expansions. 
Suppose / G Mk has g-expansion 

oo 

/(?) = J2 anqn - 

n=0 

For every prime £, there is a Hecke Operator Tg : Mk — >■ given by 

oo oo 

(2i/)(g) = ]T a n ,g™ + ^ ^ a n q nl . 

n=0 n=0 

An important map is multiplication by the Hasse invariant A, dehned in 2.1. As we 
mentioned above, A has g-expansion A(q) = 1. Multiplication by A is an injective linear 
map 

M k — ► Mfc+^i), ft-^Af. 

Of course, it behaves like the identity map on the level of g-expansions, and therefore com- 
mutes with the Hecke Operators T%. 

If / is a modular form (mod p), its filtration is dehned by 

= min{Ä; G N | / = v4*g for some 5f G M fc , z G N}. 

2.3. The algebra of modular forms. The product of a form of weight k\ and a form of 
weight ^2 is a modular form of weight k\ + ^2- We take this multiplicative structure into 
account by setting 

M = 0M fc . 

fcez 

This is a graded F p - algebra of Krull dimension 2. The g-expansion map 

M— >F p [[g]], f^f(q) 
is an algebra homomorphism with kernel — 1)M. 



2.4. The theta Operator. There is a derivation on M, raising degrees by p + 1: 

d 

c?g" 

whose effect on g-expansions is 



^ : M fc — > M fc+(p+1) , / i— )• g— /, 



(#f)(q) = J2na n q n . 



n=0 



Katz gave a geometric construction of this Operator and described some of its properties 
in [Kat77]. Of these, we will need 

Proposition 1 (Theorem (2) and Corollary (5) of [Kat77]). 

(a) If f G Mk has filtration k and p does not divide k, then df has filtration k+p+1. 



(b) If f E Mfe has = ; £/ien / /ios o unique expression of the form 

f = A r g p , 

where < r < p — 1, r + k = (mod p), g E and p£ + r(p — 1) = fc. 

We use this to find out whether an eigenform can be in the kernel of d: 

Proposition 2. If f is a Hecke eigenform and = for some i, then f is a scalar 

multiple of some power of the Hasse invariant A. 

Proof. We start by proving the case i — 1. 

By Proposition 1, the g-expansion of / E ker$ is of the form 

f(q) = clq + a P <f + a 2p q 2p + ■■■ 

Since / is an eigenvector for T p (say with eigenvalue a(p)), we have 

a(p)a + a{p)a p q p + . . . = a(p)f(q) = (T p f)(q) = a + a p q + . . . 

We conclude that a p = 0, but then a np = for all n > 1. So the g-expansion of / is actually 
constant f(q) = üq. We normalize / so that f(q) = 1. Then A — f is in the kernel of the 
g-expansion homomorphism, so 

N 

A - f = (A - l)h for some h — ^ hj E M, 

3=0 

where hj is homogeneous of degree j. 
We distinguish three possibilities: 

(a) The weight of / is p — 1. Then / and A are both in M p _i and have the same 
g-expansion, so by the g-expansion principle f = A. 

(b) The weight of / is < p— 1. Then comparing the highest degree terms in A— f = Ah — h 
we see that A = Ah^, which means that h = 1 and / = 1. 

(c) The weight of / is > p — 1. By looking at the highest degree terms in — f+A = Ah — h 
we get / = —Ah N . Note that = $(/) = $(/t/v) an d fojv is a Hecke eigenform with 
weight strictly less than the weight of /. We repeat the whole argument with / 
replaced by Hn, until we fall in one of the cases (a) or (b), and we are done since each 
step peels off a factor of —A. 

To finish the proof, we need to consider the case i > 1. So suppose ■#'(/) = 0, and let 
g = $ l ~ 1 (f). Suppose g ^ 0, then g is a Hecke eigenform satisfying 'ß(g) = 0, so by the case 
2 = 1 proved above, we know that g = cA n for some c,n. However, since i > 1, g is in the 
image of hence g = cA n is a cusp form, which implies that g = 0. We can therefore move 
all the way down to $(/) = 0, from which we conclude by using the case % = 1. □ 

2.5. Hecke eigensystems. In view of our interest in Galois representations unramified 
outside p, we define the (away-from-p) Hecke algebra by 

^ = Z[T £ \t^p\. 

By a Hecke eigensystem we will mean a ring homomorphism 

$: — >F P . 
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It is clear that the Spaces are F p J^-modules. We say that an eigensystem $ occurs in 
Mk if there exists a nonzero / G such that 

Tf = $(T)/ for all T E Jif. 

Ii $ is an eigensystem, we dehne the (hrst) twist of $ by 

$[1] : — T £ ^£$(T,). 

It is clear that this Operation can be repeated (at most) p — 1 times before getting back to <J>. 
The resulting eigensystems are called the twists of $. The twisting Operation has a modular 
Interpretation: if $ is the eigensystem of / e and $/ 7^ 0, then $[1] is the eigensystem 
of tf/e M fcfp -,. 

We will say that two eigensystems $ and ^ are equivalent (write $ ~ \E0 if $ is a twist of 
i.e. if there exists % such that $ = 

One of the crucial results for our computational work is due to Ash and Stevens (Theorems 
3.4 and 3.5 in [AS86]): 

Theorem 3 (see Theorem 3.4 in [Edi92]). Every modular eigensystem has a twist that occurs 
in weight < p + 1 . 

This indicates that, instead of having to work with spaces of arbitrary weight, it suffices 
to restrict to weight < p + 1 and take twists. 

2.6. The Sturm-Murty bound. We need to be able to decide whether two eigensystems 
are equal by comparing only finitely many of the eigenvalues. The following result (due to 
Sturm and revisited by Murty) solves this problem in the case of two eigenforms of the same 
weight: 

Theorem 4 (special case of Theorem 1 in [Mu 7]). Let f and g be holomorphic modular 
forms of weight k and level 1, with Fourier coefficients a/(n) and a g (n). Let ß{k) = k/12 
and suppose that 

üf(n) = a g (n) for all n < ß(k). 

Then f = g. 

The proof works in any characteristic; via the relation between Fourier coefficients and 
Hecke Operators we arrive at the form in which we will use the result: 

Proposition 5. Let $ and \& be eigensystems occurring in the same weight k and suppose 
that 

$(£) = \]/(£) for all primes t < ß{k). 

Then $ = ^. 

3. SOME CONSEQUENCES OF THE THEORY OF THETA CYCLES 

Let / be a modular form such that $/ 7^ 0. The -Q-cycle of / is defined to be the (p — 1)- 
tuple of integers 

A lot is known about the structure of theta cycles. For low weights, we will use the 
following Classification given by Edixhoven: 
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Proposition 6 (Proposition 3.3 in [Edi92]). Let p > 5 be prime. Let f be an eigenform 
(mod p) of weight and filtration k, where k < p+ 1. Let (ae) denote the eigenvalues of f . 
(a) If a p 7^ (f is ordinaryj, then the d-cycle of f is given by 



weight 


d-cycle 


4<k<p-l 
k = p + 1 


(k + (p+l),... 
k' + (p + l),... 
(p+l + (p+l), 


k+(p-k)(p+l), 
,k' + {k-l)(p + l)) 
...,p + l + (p-l)(p + l)) 



where k' = p + 1 — k. See Figure 1. 
(b) If a p — zs non-ordinaryj, i/ien £/ie d-cycle of f is given by 



weight 



4 < k <p- 1 
= p + 1 



d-cycle 



(k + (p + 1), . . . ,k + (p - k)(p + 1), k", 
k" + (p + 1), . . . ,k" + (k - 3)(p + l),k) 
does not occur 

where k" = p + 3 — k. See Figure 2. 

Remark. We have extracted from the statement of Proposition 3.3 in [Edi92] only the parts 
that are relevant to level 1. We have also eliminated the unnecessary requirement that / be 
a cusp form (see Section 7 in [Joc82]). 

Lemma 7. Let f\ and f'2 be eigenforms with $(/i) ~ $(/2)- Then the d-cycles of fi and $2 
are the same up to a cyclic permutation. 

Proof. We start by reducing to the case where neither f\ nor f2 is in the kernel of d. Suppose 
fi G ker(?9), then by Proposition 2 we know that fi = cA n for some c, n. Therefore $(/i) = 
<3>(t4) = $(G p+ i)[p — 2], so we may replace fi by G p+ i, which is not in the kernel of The 
same goes for f 2 . 

Since the eigensystems are equivalent, there exists an integer i such that $(/i) = «ä?^/^)- 
In particular, the weight of fi and the weight of ^$2 are congruent modulo p — 1. We have 
that $(fi) 7^ and 'Oiß 1 $2) 7^ 0, so ^(/i) and 'd l+1 {f2) have the same g-expansion, and their 
weights are congruent modulo p — 1. Without loss of generality, the weight of is less 

than or equal to the weight of q & l+1 {f2)- I so there exists j such that A^(f\) has the same 
weight as $ l+l {f2)- These forms also have the same g-expansion, so they must be equal: 

But then for all a > 1 we have 

A j $ a h = $ i+a f 2 . 

Since w(Ag) = w(g) for all modular forms g, we conclude that the $-cycles of f\ and f2 are 
the same up to a cyclic permutation. □ 

We use Edixhoven's result to determine when two eigensystems are equivalent, and to 
estimate the number of twists of a given eigensystem. 

Theorem 8. Let f\, respectively f2 be eigenforms of weight and filtration k\, respectively k2, 
where 

1 < k x < k 2 < p+ 1. 
Suppose $(/i) ~ ^{f'2), then we must be in one of the following two situations: 
(a) a p (fi) 7^0^ a p (f 2 ) and k x + k 2 = p + 1; 




•fc + (p- l)(p + l) 



FlGURE 1. Theta cycles of ordinary forms: 4 < k < p — 1 (left, k' = p + l — k), 
k = p + 1 (right). 

(b) a p (/i) = = a p (f 2 ) and k x + k 2 = p + 3. 

Proof. By Lemma 7, the ^-cycles of /i and /2 are the same up to a cyclic permutation. The 
two cases now follow by comparing the general shape and the low points of the cycles in 
Edixhoven's Classification. □ 

Proposition 9. Let f be an eigenform of weight and filtration k, where 1 < k < p + 1, and 
let $ = $(/)• Let n($) denote the number of distinct twists o/$. Then n($) — p — 1 unless 
a p 7^ and k = in which case 

n(*)e{^,p-l}; 

Proof. Suppose n($) ^ p — 1. Then n($) is a divisor of p — 1, and the -#-cycle of / consists 
of copies of subcycles of length n($). 

According to Edixhoven's Classification, the only cycle that could have this shape occurs 
when a p ^ 0, 4 < k < p — 1. This shape has two low points, so > (p — l)/2. Moreover, 
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FlGURE 2. Theta cycle of a nonordinary form; 4 < k < p — 1, k" = p + 3 — k. 
the two low points must agree, i.e. 

P+ 1 

Ä; + p + l = 2p + 2-Ä; fc = f-y- . 

□ 

An immediate consequence is: 

Corollary 10. Ifp = 1 (mod 4) zs a prime, then every Hecke eigensystem modp has exactly 
p — 1 twists. 

Example. In Section 4 we prove that if p = 3 (mod 4), G( p +i)/2 always has ^-cycle of length 
(/' D2. 

If / is a cusp form of weight (p + l)/2, its ^-cycle length can be either (p — l)/2 or p — 1. 
We give an explicit example for each of these two cases. 

(a) The smallest example of a cusp form of weight (p + l)/2 with theta cycle of length 
(p- l)/2 is A mod 23: 

A(g) = g + 22g 2 + 22g 3 + g 6 + g 8 + 22g 13 + 22g 16 + g 23 + 22g 24 + g 25 + 0(g 26 ). 

I claim that t? 12 A = A 12, dA, and hence the $-cycle of A has length 11. This alleged 
equality takes place in weight 300, where the Sturm bound is 25, so it sufhces to check 
it on g-expansions up to that precision: 

(tf 12 A)(g) = g + 21g 2 + 20g 3 + 6g 6 + 8g 8 + 10g 13 + 7g 16 + 22g 24 + 2g 25 + O(g 26 ), 
(A 12 tfA)(g) = g + 21g 2 + 20g 3 + 6g 6 + 8g 8 + 10g 13 + 7g 16 + 22g 24 + 2g 25 + O(g 26 ). 

(b) The smallest example of a cusp form of weight (p+l)/2 with theta cycle of length p— 1 
occurs for p = 43. The space of cusp forms of weight 22 is one-dimensional; denote its 
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normalized generator by A 22 (an explicit expression for it is A 22 = 41G 4 G6 + 18G , 4Gg). 
The beginning of its g-expansion is 

A2 2 (g) = q + 13g 2 + 27g 3 + 41g 4 + 39g 5 + 0(g 6 ). 

The following shows that the $-cycle length is not 21: 

(tf 22 A 22 )(g) = g + 13g 2 + 4g 3 + 18g 4 + 16g 5 + 0(g 6 ), 
(A 22 tfA 22 )(g) = g + 3g 2 + 12g 3 + 3g 4 + llg 5 + 0(g 6 ). 

4. Eigensystems Coming from Eisenstein series 

Proposition 11. Let 4<ki<k 2 <p+l and let $i, $ 2 denote the eigensystems of the 
Eisenstein series Gk ± and Gk 2 - Then $i ~ $ 2 if and only if ki + k 2 = 2 (mod p — 1). In 
i/ws case, $ 2 = $i[p — fcjj. 

Proof. Suppose fci + fc 2 = 2 (mod p — 1). On one hand we have 

$ x [p - fci](T^) = f"* 1 (1 + £ fcl - 4 ) = £ p ~ kl + 1. 

On the other hand: 

fci + fc 2 = 2 (mod p — 1) =>- fc 2 = p + 1 — fci (mod p — 1), 

so 

$ 2 (r £ ) = i + = i + e p+1 - ki -\ 

For the other implication, suppose $ 2 = for some z. This means that 

f + t +kl ~ l = 1 + (mod p) 

for all primes £ ^ p. Let a, b, c be the respective remainders of the division by p — 1 of i, 
% + ki — 1, k 2 — 1. (In particular, a,b,c < p — 1.) Then in F p we have 

(3) a a + a b = l + a c for all a G F p x . 

Consider the polynomial 

f(x) =x a + x b -l-x c e¥ p [x}. 

The degree of / is at most p — 2 (or / is the zero polynomial). If / ^ then / has at most 
p — 2 roots in F p . However, equation (3) implies that / has p — 1 roots in F p , so we must 
have that / = 0. 

We have two possibilities: (i) a = and b = c, which implies % = and k\ = fc 2 , 
contradicting the assumption that k\ < k 2 ; (ii) 6 = and a = c, which implies 

k\ + k 2 = 2 (mod p — 1) and i = k 2 — 1 = p + k 2 — 2 = p — fei (mod p — 1). 

□ 

Proposition 12. Let 4 < k < p + 1. The Eisenstein series Gk has p — 1 twists, unless p = 3 
(mod 4) and = (p + l)/2, m which case Gk has (p — l)/2 twists. 

9 



Proof. The case p = 1 (mod 4) is taken care of by Corollary 10. 

Suppose now that p = 3 (mod 4). Eisenstein series are always ordinary, so a p ^ 0. If 
k (p + l)/2 then £^ has theta cycle of length p — 1 by Proposition 9. In the remaining 
case k = (p + l)/2, let $ be the eigensystem of G^. We easily see that 

$(2>) = i + £(p+D/2-i = ! + £( P -D/2 

$[(p-l)/2](T € ) = £(P-l)/2 ( 1+ £(p-1)/2) =£ (p-l)/2 + 1; 

so $ has (p — l)/2 twists. □ 

Corollary 13. The number of distinct eigensystems (mod p) coming from Eisenstein series 
is (p - l) 2 / 4 - 

Proof. This follows via simple arithmetic from Propositions 11 and 12. □ 

We end this section by discussing the possibility that an Eisenstein series and a cuspidal 
eigenform of small weights have equivalent eigensystems: 

Proposition 14. Let Gk be the Eisenstein series of weight k < p + 1 and fix an even integer 
k' 7^ 14 with 12 < k' <p+ 1. A cuspidal eigenform f of weight k' with Q(Gk) ~ &(f) exists 
if and only if k' = k and p divides the numerator of the k-th Bernoulli number Bk- 

Proof. The argument can be extracted from page 334 of [Scr73]; we include it here for 
completeness. 

Suppose there exists a form / with the given properties. Then there is some integer i 
such that $(/) = $(Gfc)[i], i.e. $/ = d l+l Gk- The conditions imposed on k' exclude the 
possibility of it being divisible by p, therefore the filtration of $/ is k! +p + 1. Similarly, the 
filtration of d l+l Gk is k + (i + l)(p + 1). Therefore 

k' + p + l = k + (i + l)(p + 1). 

However, k'<p+lsok' + p + l< 2(p + 1), from which we conclude that i = 0, so k' = k. 

Therefore d(f — Gk) = 0. Again since k is not divisible by p we get that / = G^, in 
particular the constant term of Gk is zero; but this constant term is the reduction modulo 
p of Bk/(2k), therefore p must divide the numerator of Bk/(2k). Using one last time the 
condition k < p + 1 we conclude that p divides the numerator of Bk/(2k) if and only if it 
divides the numerator of Bk- □ 

5. BOUNDS ON THE NUMBER OF EIGENSYSTEMS 

In this section, we derive an explicit formula for the well-known upper bound on the 
number 2 N(2,p) of level 1 Hecke eigensystems modulo p. 

Let iVt„ist(2,p) be the number of equivalence classes up to twist of level 1 Hecke eigensys- 
tems modulo p. We have seen that any eigensystem has at most p — 1 twists, so we get the 
inequality 

N(2,p)<N twist (2,p)-(p-l). 

2 We use Khare's notation, which is motivated by the fact that this is the number of continuous semisimple 
odd representations 

p: Gal(Q/Q)^GL 2 (F p ) 

that are unramified outside p. Note that we do not restrict our attention to irreducible representations here, 
but by Corollary 13 the difference is known to be (p — l) 2 /4. 
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if p = 


1 


(mod 12) 


if p = 


5 


(mod 12) 


if p = 


7 


(mod 12) 


if p = 


11 


(mod 12) 



We know that each eigensystem occurs, up to twist, in weights at most p+1. Therefore we 
can bound iV twist (2,p) by the sum of the dimensions of the Spaces M fc (SL 2 (Z)) for k < p + 1: 

P+i 

AW2,p) < X> mM *( SL 2^))- 

k=A 

We now use the classical dimension formulas (see, e.g. Corollary 1 in Section 1.3 of [ZagOS]): 

{0 if k < or k is odd 

L^J if k = 2 (mod 12) 
|_y|J + 1 otherwise. 

After a straightforward calculation, we obtain the following expression for the sum of 
dimensions (write Q for the quotient of the integer division of p + 1 by 12): 

'3Q 2 + 4Q 

^dimM fc (SL 2 (Z)) = 

fc=4 

[3Q 2 + 3Q 

It remains to multiply this value by p — 1 in order to obtain the desired upper bound on 
N(2,p). Note that this upper bound is asymptotic to p 3 /48 as p — > oo. 

6. Detailed description of the algorithm 

Step 1: Get the eigensystems Coming from Eisenstein series. According to Propo- 
sition 11, the complete list of such eigensystems up to twist is: Gk for 4 < k < {p + l)/2, 
together with G p+ i. 

Step 2: Get the eigensystems Coming from cusp forms of weight up to p+ 1. Fix 

a weight 12 < k < p + 1. We took two different approaches: 

• Compute the (cuspidal) Victor Miller basis over ¥ p of weight k up to and including 
the p-th coefhcient, then decompose the span of this basis into Hecke eigensystems. 

• Compute the (cuspidal) modular symbols of weight k and sign +1 over F p , then 
decompose into Hecke eigenspaces. 

This gives us a list of cuspidal eigenforms fi, ■ ■ ■ , f n with n < dim^. 

Step 3: Remove duplicates (up to twist). Three special circumstances can arise: 

(a) If p divides the numerator of then one of the f/s is the Eisenstein series Gk, by 
Proposition 14. Determine the relevant fj (need to compare g-expansion up to ß{k)) 
and remove it from the list. 

(b) If k > (p + l)/2: for each j G {1, . . . ,n} such that a p (fj) ^ 0, there could be a 
companion form in weight p + 1 — k. In this case the eigensystem of fj is a twist of 
an eigensystem that has already been listed, so we should remove fj from the list. 
Checking this requires comparing the ordinary /, 's with the ordinary forms of weight 
p + 1 — k, up to precision ß{k +p + 1). 

Here is the justification for the comparison bound: we have / of weight k > (p+l)/2 
and g of weight p + 1 — k. We want to check whether the g-expansions *&f (in weight 
k + p + 1) and d k g (in weight kp + p + 1) are equal. A priori it seems that this must 
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be checked in weight kp + p + 1, where we are verifying the equality A k $f = d k g. 
However, as Buzzard pointed out to us, we can do much better by using theta cycles. 
We are in the Situation illustrated in Figure 1: i?/ is the first low point of the cycle, 
and 6g is the second low point. Following the cycle, we see that $ h g is back at the 
first low point, i.e. that d k g has filtration k + p + 1. Therefore it suffices to perform 
the comparison in weight k + p + 1, checking g-expansions up to ß(k + p + 1). 
(c) If k > (p + 3)/2: for each j G {1, . . . ,n} such that a p (fj) = 0, there exists a non- 
ordinary form g of weight p + 3 — k with the same eigensystem up to twist, therefore 
fj should be removed from the list. 

As a consistency check (inexpensive since there are not many nonordinary forms), 
we could actually compare the nonordinary fj's to the nonordinary forms of weight 
p + 3 — k and find the corresponding form there. As we see from the theta cycle of 
/ in Figure 2, we want to check whether the g-expansions of d p+l ~ k f and g agree, in 
weight p + 3 — k. We therefore need to compare up to precision ß(p + 3 — k). 
We now have the list of all eigensystems up to twist. 

7. SUMMARY AND DISCUSSION OF RESULTS 

The table appearing below in the appendix records, for all the primes under 2000, the 
number of distinct non-Eisenstein 3 eigensystems mod p, the upper bound on this number, 
as well as any interesting features that each prime might have. The latter are denoted by 
an E for Eisenstein-cuspidal congruences, a C for companion forms, or an N for nonordinary 
forms, followed by the weights in which the corresponding phenomenon occurs. Note that 
companion forms and nonordinary forms always show up in pairs, but only the smallest 
weight is recorded in the table for each such pair. 

The first explicit examples of companion forms appear in [Gro90], resulting from compu- 
tations done by Elkies and Atkin. They focused on finding primes at which the reduction 
of the six cuspidal eigenforms with rational coefficients have companions. Higher degree 
examples were given by Centeleghe in his thesis [Ccn09], going up to p = 619. Our results 
extend this ränge to all p < 2000. 

Similarly, we find new examples of nonordinary forms mod p < 2000 of weight k < p + 1, 
extending those listed in Tables 5 and 6 of [Ccn09] and the results of Gouvea in [Gou97] . 

It is interesting to compare our results with Centeleghe's table in [CenlO]. Out of the 299 
lower bounds he computes, 164 are marked with a star, meaning that they are proved to give 
the actual number of representations. Our results indicate that a further 111 of his lower 
bounds coincide with the exact numbers, for a total of 275 out of 299. 

Finally, we notice that the "interesting" phenomena described above are quite rare, and 
the actual number of eigensystems deviates very little from the explicit upper bound given 
in Section 5. For instance, among the last 20 primes in the ränge we computed, the relative 
difference between the actual number and the upper bound is always less than 0.017%. This 
relative "error" is plotted in Figure 3 at three different zoom levels. Note that the primes 
congruent to 1 mod 4, represented by blue discs in the figure, generally tend to be closer to 
the upper bound than the primes congruent to 3 mod 4. 



We decided to exclude the Eisenstein eigensystems from the count in order to ease comparison with 
Centeleghe's results. As Corollary 13 indicates, the number of Eisenstein eigensystems (modp) is (p— l) 2 /4. 

12 




800 1000 1200 1400 1600 1800 2000 



1000 1200 1400 1600 1800 2000 



FlGURE 3. The relative difference (as a percentage) between the actual num- 
ber of eigensystems and the upper bound. The three different views show 
the primes: (top) 11 < p < 2000; (middle) 500 < p < 2000; (bottom) 
1000 < p < 2000. The blue discs represent primes congruent to 1 mod 4, 
while the red triangles represent primes congruent to 3 mod 4. 
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